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1 Introduction 

It is well-known that the geometric description of the Euler-Lagrange equations of a me- 
chanical system, derived from Hamilton's principle, heavily relies on the intrinsic geometry 
of the tangent bundle TQ of the configuration space Q. In the case that the Lagrangian 
L is invariant under the action of a Lie group G, Hamilton's principle can be reformulated 
into a reduced variational principle that takes into account the symmetry properties of 
the system. The advantage is that the reduced variational principle leads to equations 
defined on a reduced space, i.e. equations depending on 'fewer independent variables'. 
Thus, instead of considering the Euler-Lagrange equations on the total space Q, we are 
interested in the reduced equations, the so-called Lagrange-Poincare^ equations, which 
are formulated on the Atiyah quotient bundle TQ/G Q/G (see e.g. [15 )• It is clear 
that the geometry of such quotient bundles now becomes of interest. Weinstein has 
pushed ahead our understanding in this matter by showing that the geometric structure 
which lies at the heart of the Lagrange-Poincare equations is essentially the same as the 
one of the Euler-Lagrange equations, namely that of a Lie algebroid. Therefore, the ge- 
ometry of 'a Lagrangian system on a Lie algebroid' unifies the geometry of both standard 
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Lagrangian systems and those where the symmetry properties are taking into account. In 
the case of standard Lagrangian systems, the tangent bundle carries a canonical Lie alge- 
broid structure which is given by the usual Lie algebra of vector fields on Q. For systems 
with symmetry, we make use of the so-called Atiyah algebroid structure on TQ /G Q/G 
to describe the evolution equations (see e.g. [Ul ITU]). 

In this paper, we mainly deal with mechanical systems that are subject to some kinemat- 
ical constraints (i.e. depending on the velocity), also called non-holonomic constraints. 
The governing equations for these models are the so-called Lagrange-d'Alembert equa- 
tions. For systems having additional symmetry properties, the equations of motion can 
be reduced to the Lagrange-d'Alembert-Poincare equations (a fairly complete list of ref- 
erences on non-holonomic mechanics can be found in PI IE])- The main purpose of this 
paper is to show that the Lie algebroid structure constitutes a unifying geometric struc- 
ture for describing simultaneously both types of Lagrangian systems with non-holonomic 
constraints. Throughout this paper, we will develop the basic geometric concepts and 
objects that are involved. 

We believe that at this stage it is instructive to provide a local version of what will follow 
(we assume that the reader is familiar with the concept of a Lie algebroid). We recall the 
equations of motion for a non-holonomic mechanical system and take them as the starting 
point for further generalisation in the framework of Lie algebroids, eventually leading 
to the equations describing a 'non-holonomic Lagrangian system on a Lie algebroid'. 
Let L G G°°{TQ) be the Lagrangian of a non-holonomic system. The constraints are 
assumed to define a subbundle of the tangent bundle: t : D ^ TQ, i.e. locally we have 
that l\{x)w^ = v'' where v'^ denote the components of i{w) w.r.t. the standard basis 
{d/dx''} and where are the components of w w.r.t. some basis of D. A basis of 
the annihilator space is denoted by = uj\dx'^, i.e. oj\i\ = 0, \/b,A. The Lagrange- 
d'Alembert principle states that the equations of motion for such a non-holonomic system 
are determined by 5 J L{x^{t),x^{t))dt = 0, where the variations 5x* should satisfy the 
constraint, i.e. 6x G -Da,(t) for each t G [a,b], and moreover 6x{a) = 5x{h) = 0. The 
induced equations are called the Lagrange-d'Alembert equations: 



for some functions Xb{t). It is not hard to see that this system is equivalent to the system: 



Next, we repeat the above construction for systems on a Lie algebroid. Assume that a local 
coordinate chart {U, (x, v)) of a Lie algebroid r : V ^ M is given. Given a Lagrangian L 




(1) 
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on V, the Lagrangian equations on the Lie algebroid are given by: 

X* = pjj(x)v", 

d_fdL\ ^^dL_,dL (3) 
dt\dv-J ^"^x* dv^' 

where and C^^, are the structure functions of the Lie algebroid. We now assume that 
W is a subbundle of V with injection i : W ^ V. We will use the observation of the 
previous paragraph: by contracting the unconstrained Lagrangian equations with the 
components of the injection l'\ of the constraint distribution into the tangent bundle, 
the non-holonomic equations can be produced. Assuming that Lagrangian systems con- 
strained to the subbundle W of a Lie algebroid V have a similar behaviour, we can now 
postulate that the constrained Lagrangian equations are: 

Pl{x{t))y%t), 

We can now write the above equations in terms of derivatives of the constrained La- 
grangian, defined by Ld^x^, w"^) = v° = i'\w"^). Then, along the solution v = iow of 
(0)), a straightforward calculation leads to: 

d_ fdLA 



x'{t) 
v'^(t) 



d ( dL 
di [dy^ 



_ d dL d I dL\ 

~ dx^ dv-^^y^'^dx^ dv-j 



If we define = p'J'X then, since dLc/dx"^ = dL/dx'' + {dL/ dv'^){di% / dx'^)w^ , it is clear 
that the curve {x(t),w(t)) is a solution to the system: 

{x* = X\{x)w^, 
_ / 5^ (5) 

dt [dw^J ~ dx^ + [<-bJB^A Q^, + ) ■ 

In the first part of this paper, we develop all geometric structures required to provide an 
intrinsic formulation for this system of equations in (jSJ. In ^3], E. Martinez presented a 
solid geometrical framework for Weinstein's systems on Lie algebroids. His approach is 
very similar to the usual formalism for Euler-Lagrange equations. The important differ- 
ence is, however, that no longer a vector field is the main geometrical object, but rather 
a section of a 'prolongation bundle' inducing a unique vector field on the Lie algebroid. 
In this paper we extend the framework of Martinez to the above constrained systems 
(0). For that purpose, we first define in Section |21 an exterior derivative on a subbundle 
of a Lie algebroid. Next, we develop the concept of prolongation bundles in Section 01 
eventually leading to all necessary tools for an intrinsic description of (0) in Section EJ In 
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our formalism, the system (0) will be regarded as a section of an appropriate prolongation 
bundle. With this section a vector field is associated, whose integral curves are precisely 
the solution of (0). The second part of the paper is devoted to examples from known dy- 
namical systems that allow a formulation in terms of the above equations (cf. Section [^I). 
These examples only deal with autonomous systems with linear constraints. Finally, we 
will discuss some of the advantages of our approach and we will indicate some directions 
for future work. 



2 Exterior derivatives 

A Lie algebroid is a vector bundle t : \J ^ M with a real Lie algebra bracket on its 
set of sections [■, ■] : Sec(r) x Sec(r) Sec(r). Moreover there is a linear bundle map 
p : y TM (and its natural extension p : Sec(r) —>■ X{M)) which is related to the 
bracket in such a way that, for all s, r G Sec(r), / G C°°{M) 

[s,/r] = /[s, r]+p(s)(/)r 

is satisfied. 

Choose a basis {ea} for Sec(r) and denote the corresponding local coordinates on V by 
(a;*,v"). Then, the structure functions of the Lie algebroid are smooth functions p^ and 
C^j, on M which satisfy p{ea) = Pa-£j and [ea, e^] = C^ah'^c- In Lie algebroid theory, the 
role of differential forms is played by sections of exterior powers of the dual vector bundle, 
i.e. skew-symmetric, C°°(M)-linear maps lo : Sec(r) x ■■■ x Sec(r) — > C°°(M) (with h 
arguments) will be called k-forms on Sec(r) and the set of all such forms will be denoted 
by /\'^{t). The defining properties of a Lie algebroid structure lead to the definition of an 
exterior derivative on /\{t). Let cu G /\'^{t), then the (k + l)-form dcu is given by 

fc+i 

du;(si, . . . ,Sfc+i) = ^(-l)*"V(si)(^^(si, • • • ,Si, . . . ,SA,.+i)j 

i=l 

+ Yl (-l)*^^^([si,Sj],Si,...,Si,...,s},...,Sfc+i). (6) 

i<i<j<k+l 

The operator d has the property d{u!i A UJ2) = dui A 002 + {—l)^'^uji A duj2, with ki the 
degree of ui, and moreover d"^ = 0. Locally, dx^ = pj^e" and de'^ = — |C^b^° ^ where 
the set {e"} stands for the basis of Sec(r*) which is dual to {ea}- More details on the 
properties of the exterior derivative can be found in ^Hl E] • 

Suppose now that a vector subbundle /i : W ^ M of r is given. It is obvious that each 
coordinate system {x\ w^) on W can be extended to a coordinate system (x*, v"^ = w^, v"") 
on V. However, in what follows, it will be more convenient to consider the basis {6^} of 
Sec(r) to be a priori given. Then, the sections of an arbitrary basis {e^} of Sec(p), can be 
written as = i^^a injection can be denoted by i : W — ^ V; (x*, w^) ^ {x\ v° = 

i'\u)'^). We will use A for the restriction of p to W: A = p o i ; W — > TM; {x\w^) t— > 
{x\ X* = A^w^ = pl^i'Xw^). 
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Forms on Sec(r) can be pulled back to forms on Sec(/i). Indeed, if is a fc-form then i*uj, 
defined by 

t*uiWi, ...,Wk) = toim,, ...,iWk), WiG Sec(/x) 

is a fc-form on Sec(r). By composing the exterior derivative d with i* we can define a 
mapping from /\''{t) to A'^^^I/^)? which is denoted by S. Thus, if is a fc-form on Sec(r), 
then 

6u = i*duj, (7) 
is a (A; + l)-form on Sec(/i). The operator 6 satisfies the rule 

6{uJi A UJ2) = 6uJi A i*uj2 + {-l)''H*uJi A 6002, 

and could be called a derivative along i for this reason, but we will simply refer to 6 as 
an exterior derivative. Obviously, 6 o d = 0. In the above introduced coordinates, 

6x' = pii'Xe^ = X\e^ and 5e^ = -^D^^e^ A e^, 
where 0%^ = C'^Aic- 

To illustrate the above notions, we look at the case where W is a Lie subalgebroid of V. 
A more general definition of a Lie subalgebroid W — of V — M can be found in ^T] . 
Here, we will only consider the case that the base manifolds coincide. A Lie algebroid 
/i : W ^ M (with anchor map A and structure functions -D^^) is a Lie subalgebroid of 
r : V ^ M (with anchor map p) if there exist a injective morphism i : \N —>■ V such that 
X = p o i and 

i[W,, W2] = [liWi], tiW2)], e Sec(/x). (8) 

In coordinates the above expression reads: 

>cfSB = i'.«-Ai§ + A^§. (9) 

There are now three exterior derivatives around: d, 6 and the exterior derivative d^ on 
the Lie algebroid p. The condition (jH)) can equivalently be written as 

(l^^oi* = i* od = 6. (10) 

Indeed, for functions on M it is obviously satisfied, while for 1-forms 6 G /\^{t) 

d^'eeiWu W2) = x{Wi)e{iW2) - x{W2)e{Wi) - e{i[Wu W2\). 

This is d6{iWi,iW2) if and only if (jSJ is satisfied. The proof then follows from induction. 
It is now also obvious that for Lie subalgebroids 

d'' 06 = i* od'^ = {d^ f oi* = 0. 
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3 Prolongation bundles 



As stated in the introduction, an important role in our description will be played by the 
so-called prolongation bundle associated to a Lie algebroid. Sections of such bundles will 
eventually lead to the vector field on the Lie algebroid that generates the equations in © . 
Furthermore, these bundles allow a generalisation of all intrinsic objects defined on the 
tangent bundle needed to write down the standard Lagrange equations. We start with 
defining the prolongation of an arbitrary bundle, and then continue with the introduction 
of all the structures that lead to an intrinsic formulation of (0). 

Let IT : E M he an arbitrary fibre bundle. The prolongation of tt by an anchored 
vector bundle r : V ^ M is the (vector) bundle irf : T^E E. Here, the total space 
TPE is the pullback manifold p*TE = {(v,Xe) G V x | p(v) = T7r(Xe)}. The natural 
projections of T^E onto V and TE will be denoted by vr^ and p'^, respectively (see the 
diagram below). On the other hand, the bundle projection ir^ is given by Te o p^, i.e. 
7Tp{v,Xe) = e (see also e.g. HDl El CHI EI). 




Suppose we fix a bundle adapted coordinate chart (f/, {x\y°')) about a point e E E. It is 
not difficult to see that {Xa,ya}, with 



d 

Xa{e) = { e,(7r(e)),p^(x) — 



and 



0(vr(e)), 



d 



is a basis of Sec(7r''), induced by the basis {ea} of Sec(r) and the basis ^} of X{E). 
The set V^E = {(0, X^) G T^E} is called the bundle of vertical elements of T^E. A basis 
for the set of vertical sections Ver(7r'') is given by the {D^o}. In what follows, we will use 
{A"^, 3^"} for the basis of Sec((7r'')*) that is dual to {Xa, 3^a}. 

Suppose now that r : V — M is a Lie algebroid. The Lie algebroid structure on r can be 
naturally extended to a Lie algebroid structure on the prolongation bundle ttP : T^E E. 
The anchor map of this Lie algebroid is p'^ : T^E TE, with 



Pa 



d_ 



and 



_d_ 
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and the bracket is given by 



[^a. Xb] = ci,x,, \x^, yp\ = 0, [y^, y^] = o 

(for more details, see e.g. [El USD- 
Let n : \N —y M he Si subbundle of r. We will need three different prolongation bundles. 
The first is the one where tt is also r : V M, that is : T^'V — * V. Since this bundle 
has a Lie algebroid structure, it also has an exterior derivative d, locally characterized by 
the relations 

dx' = p^X", dy" = V^ dX' = -^C'^f^X^AX^ and dV = 0. 




The second prolongation of interest is fi^ : T'^W W. Again, fi^ is a Lie algebroid with 
a corresponding exterior derivative, denoted by d. If {Xa, Wa} stands for the basis (fTT| 
in this situation, then 

dx' = piX'', dw^ = W^, dX" = -^^Clf^X" A X^ and dW^ = 0. 

/i'' is in fact a Lie subalgebroid of r'' (in the more general sense of [Hj, see also ^^). 

Last but not least, we will also need the prolongation p^ : T^W —>■ W, which is a vector 
subbundle of the Lie algebroid T^W — > W. Let {Xa, Wa} stand for the basis of Sec{p^), 
induced by the bases {e^} and g^}- Then, the injection is given by the map 
I : T^W TP\N; {wi, X^J ^ {i{wi), xjj which, locally, is of the form 

I{Xa) = i\Xa and /(Wa) = Wa- 

In view of the discussion of the previous section, we can use I to introduce an 'exterior 
derivative' 5, defined as in ((7j), which will map forms on Sec(/i'') to forms on Sec(/i^). 
Locally, 

5x' = \\X^, 5w^ = W^, SX' = -^D'XbX^AX^ and SW^ = 0. 

7 



A second important map T^i : T^'W — > T^M is the linear bundle map over i : W ^ V, 
given by T^iiy^Xyj) = (v, Tj(^)i(X^)). Let an element Z G Z^W be given by Z = 
X'^Xaiw) + W^Wa{w), then 

Both prolongations : T^W V and jj^ : T^W ^ W carry further interesting canonical 
objects. For example, there exists a naturally defined vertical lift ^ : r*V —>■ VV C T^V. 
Indeed, for (a,v) G r*V, we can define first a vertical element G TgV by means of its 
action on functions / G C°°(V), 

Then, the element v^' G VV is defined as (0, v^). We have = Va- Evidently, there 
exists a similar notion for the prolongation /x'**, here the vertical lift is a map ^ : /i*W ^ 
V'^W C T^W. Next, there are also two so-called vertical endomorphisms S'^ = ^ o : 
Sec(r'') Ver{Tp) and = ^ o i2 '■ Sec(/i'^) — > Ver^p}), where ji stands for the 
projection (v2,X^J G T^'V t-^ (vi,V2) and ja : /i*W can be defined analogously. 

In the above introduced bases, 5"^ = X"- ® Va and = ® W^. Remark that, 
although the projection : T'^W — > /i*V can easily be defined, there exists no vertical 
endomorphism on T^'W because there is no vertical lift which maps elements in /i*W to 
VW. Two final important objects are the Liouville sections = v'^Va G Sec(r'') and 
C = w^Wa e Sec(/i^). 



4 Lagrange equations on a subbundle of a Lie alge- 
broid 

In this section we formulate the equations for a non-holonomic system on a Lie algebroid. 
The main purpose is to arrive at an equation that uniquely determines a vector field on 
W whose integral curves are precisely the solutions to the constrained equations (0). We 
will only deal with regular Lagrangians L G C°°(V), i.e. we will assume that the matrix 
{d"^ L / dy°' Om^) is regular at every point. 

Let's recall first briefiy Martinez's description of Lagrangian systems on a Lie alge- 
broid. An important subclass of sections of the prolongation bundle r'' : T^'V V are 
the so-called jJsendo-SoDEs jTH HH], or simply SoDEs in [HIl US]- They are sections V 
of such that r2 or = id (r^ is the projection (v, y) G T^'V v G V). Locally, a 
pseudo-SoDE is of the form 

r = y'^Xa + rVa. 

It is not difficult to see that pseudo-SoDEs are in a one-to-one correspondence with vector 
fields X on V with the property that p(v) = Tr(X(v)). Keeping this in mind, by defining 
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on V and a 1-form 



a pseudo-SoDE, one can give an intrinsic description of the equations 0- A regular 
Lagrangian L defines a function 

BT 

EL = p^{C^)L-L = y- — {y)-L{y) 

r)T 

e, = s^idL) = ^Mx^ (12) 

on Sec(r''). The dynamics are then given by a pseudo-SoDE V of the prolongation 
that solves the equation 

i^dOL = -dEL. (13) 

Solutions of the equations Q are then nothing but integral curves of the associated vector 
field p"(r) G x{y). 

For the constrained systems we wish to preserve, as much as possible, the structure 
of the equation (|T3|) . By analogy, it is clear that we should represent the dynamics by a 
section of the prolongation bundle : T^\N W. It is easy to see that pseudo-SoDEs 
r on this bundle have locally the form 

r = w^'Xa + f^VA. 

Crucial in ()13p is the exterior derivative d of the prolongation tP. Unfortunately, the 
prolongation bundle does not carry a Lie algebroid structure and therefore there is 
no available exterior derivative. The next best thing is the above introduced operator 
5. As a consequence, the analogues Ej^ and 6l of the function E^ and the 1-form 6l to 
the constrained case should be sought, respectively, among the functions on W and the 
1-forms on Sec(yU^). For El we can simply take the restriction of El to W. If we define 
the constrained Lagrangian Lc G C°°(W) as the restriction of L to W, Lc{w) = L{i{w)), 
then it is easy to see that El can also be given by 

El = yiC^iL,) -L, = nJ^'-^iw) - L,{w) G C-(W). 

The construction of 9l in (fT^ can, however, not directly be translated to fj,'^. Indeed, 
although dLc is a well-defined 1-form on Sec(/i''), there exist no vertical endomorphism 
on Sec(yu''). On the other hand, we can also start from 6Lc, which is a 1-form on Sec(/i'^), 
but its image under the vertical endomorphism S'^ of /i'^ gives us a 1-form 

Ol^ = S^iSL^) = Q^HX^ = ^\^M^))x^ 

on Sec(yU'^) and not on Sec(/i''), as required. Having found no direct construction on it 
seems appropriate to take one more step backwards (w.r.t. the diagrams of the previous 
section). We will use a suitable restriction of 6l to yU^. 

Definition 1. The Poincare-Cartan 1-form 6l is defined as the 1-form {T''i)*{9l) on 
Sec{fiP), i.e. 

eL{w){M, W) = eL{i{w)){T%y, W)), 

dL 

or, locally, Ol = ^~;7(^(^))'^"- energy function El is the restriction of El to W. 
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From the coordinate expressions it is easy to see that = I*{6l), i.e. 

eLAwi){w2,w) = eL{wi){I{w2,W)). 
We now have all the ingredients for a coordinate free description of (0). 

Definition 2. A Lagrangian system on fi^ : T^\N ^ \N is a pseudo-SODE T G Sec{^^) 
that solves the equation 

trS9L = -6El. (14) 

It can easily be checked that (|T^ gives indeed the correct equations (0). First, one can 
calculate that 



and 

_ ( ..,c-a ^-^^ I „.,B„.,c-b (^^^ ^-^^ \« _L „A-"v& 

V dx'dy'' ^ dx' dvdy^ dx' J ^ ^ " dM^'dv^ 

The coefficients /"^ of F should therefore satisfy 



The dynamics are given by the equations of the integral curves of X = A^(F) = w'^A^^ + 



f^j^, i.e. they are solutions of 



X* = X\{x)w^, 

Along such solutions, we find that 

^2 r Pit a -a . a2 ; 



(16) 



and 



(9L dL di^A . 



After plugging this information into (fT3j) . ([1^) becomes exactly (jSJ 



X* = A^(x)w^, 
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5 Examples 



The main example of this section shows that nonholonomic mechanical systems with 
symmetry admit a formulation in the above framework after reduction. The first two 
examples are straightforward. 

1. Lagrangian systems on Lie algebroids. If yU : W — > M is exactly the Lie algebroid 
r : V — i> M, then = 6^ and Lc = L. The equations (fTTj) are then exactly Q- Examples 
of such systems can be found e.g. in rigid body dynamics (see e.g. jT^I for a worked out 
example about the heavy top), control theory ^Hj or Lagrangian systems with symmetry 
on principal fibre bundles [HI • 

2. Lagrangian equations on Lie subalgebroids. With the help of 0, it is easy 

to see that in this case, (-Dr^ — A*i— + Art-^ It- — = D'^a——^. The equations (fTTIl 

V ax* ax* / o\j'^ ow^ 

then become exactly the equations for the constrained Lagrangian Lc on the Lie 

algebroid /i : W — ^ M. This is also clear from the expression ()14|1 . Indeed, in this case 

. j^Ay^ _^ yy inherits the Lie algebroid structure form /i and has also an exterior 

derivative which acts on functions and forms in the following way 

d^'x' = X\X^, d^w^ = W^, d^'X^ = ~D%X^AX^ and = 0. 

Obviously 5El^ = d^^Ei^. Moreover, due to ©, 

Therefore (jl4|) is indeed ird^Oi^ = —dEi^. 

3. Non-holonomic systems. In the introduction we already have defined the Lagrange- 
d'Alembert equations describing a non-holonomic mechanical system. We will use similar 
notations. Suppose Q is the configuration space of a mechanical system that is subject 
to some kinematic (linear) constraints. If the constraint distribution is m-dimensional, 
then, due to the regularity condition of D, we can express m velocities defined by 
X* = {r\ s°), up to a renumbering, in terms of the others 

i° = -A^(a;)r^ 



Let L : D ^ TQ, then it is given by {s°',r^,r^) ^-^ {s°',r\f\s'^ = —Af-r^). Before 
continuing, we wish to make the following remark: in some geometrical models that 
treat non-holonomic systems it is sometimes further assumed that there exists a bundle 
structure of the configuration space Q over some manifold N such that D is the horizontal 
distribution HM C TQ of a connection on Q N (this can always be done locally). 
This additional assumption is not necessary for our purposes. The purpose of this section 
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is to show that the equations for non-holonomic mechanical systems that can be found in 
e.g. P, fit into the framework presented above. 

We can take V to be TQ, equipped with the natural Lie algebroid structure: the anchor 
map p is the identity and the Lie algebroid bracket is the usual Lie bracket of vector fields. 
Then is nothing but TTQ TQ. The subbundle W is precisely the distribution 
i.e. i = L. Likewise, /i'' becomes simply TW W. Since X = i, the last prolongation 
bundle /x* (on which the dynamics will be defined!) is T*W = W Xj TW W. The 
Poincare-Cartan 1-form 6' j;^ is in this case, the puUback of the usual Poincare-Cartan form 
dL dL 

9i = ——fdr^ + 7- — ds°' for L by the injection i, i.e. it is a a differential form on W 

or^ as" 

which formally looks similar, but where the coefficients should be evaluated along the 
constraints. The biggest difference between our approach and many others is that the 
fundamental form 66 l and the dynamics F are not a differential form or a vector field on 
W, but a form and a section of a prolongation bundle. The equations (fT7|) (with = Sj 
and A} = —Aj) are of course the required Lagrange-d'Alembert equations 

d_ (^L^^ ^dL^_ dL^ _rJB-— ^^^^ 
dt\dr^ ) dr' ^^ds'^' 

dA^f dA'} .OA"} .dA^f ^, , , , , 

where Br j = -—- — — -- + ^7 -— — Aj——-. The standard way to obtam the equations 

(IIHl) goes by taking (j21) as starting point, and following the same procedure that has lead 
us to (0) in the introduction (see also equations (5.2.7) in T]). 

4. The reduction of non-holonomic systems with symmetry. In this example, we 
will call the configuration space of the mechanical system Q and the Lagrangian /. We 
further assume that / is invariant under the (tangent lift of the) action of a Lie group 
G and that tiq : Q ^ Q/G has the structure of a principle fibre bundle. We will first 
introduce the Lie algebroid structure of interest. It is assumed that the reader is familiar 
with the natural constructions associated with a principal fibre bundle |12j . 

Let Q = {Q X q) /G ^ Q /G he the associated Lie algebra bundle (for a definition see e.g. 
[T^ . it is a Lie algebroid structure with vanishing anchor map and structure functions 
the structure constants of the Lie algebra q). Suppose that A : TQ — > g is a principal 
connection on ttq with horizontal distribution H. Then TQ = H ® Vttc, where Vtiq = 
{iV'^G)q = Tq{TT~^{[q])) \ q G Q} is the vertical distribution. It is shown e.g. in 0113 that 
the connection induces a vector bundle isomorphism between TQ/G and T{Q/G) ©g 
by means of 

oiA{[vq]) = Tndvq) © [q ■ A{vg)], Vg e TgQ. 

Here, [q ■ ^] stands for the equivalence class of (g, ^) G Q x g, i.e. (g, ~ il9^ Adg-i ■ ^) for 
all g E G. Remark that aA{H/G) = T{Q/G) and aAiVira/G) = g. In the notations from 
the previous sections, M = Q/G. The Lie algebroid is then defined by r : V = TQ/G ^ 
T{Q/G) © g Q/G, where the bracket is taken to be the Lie bracket restricted to right 
invariant vector fields on Q. The anchor map p : V TM of the Lie algebroid is nothing 
but the projection onto T{Q/G). The Lie algebroid bracket can be given by 

[Xi © si, X2 © S2] = [Xi, X2] © (VxiS2 - Vx.si - cu(Xi, X2) + [si, S2]) (19) 
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with Xi G A'(M), Si G Sec(g). V denotes the covariant derivative on g — M associated 
with the connection A and uo is its curvature. 

Recall that a local trivialisation of the principal fibre bundle induces a bundle adapted 
coordinate chart on every associated bundle of ttg- In particular, if is a basis for the Lie 
algebra g, then a basis for the set of local sections of g are defined by ea(x) = e)-ea\ 
where xeUcM = Q/G and ip : 7!'q^{U) ^ U x G is such a local trivialisation of 
ttg- Suppose are coordinates on g w.r.t. this basis. Then, V has coordinates 

(a;*,v" = Using this coordinate system, the coefficients of the Lie bracket [si,S2] 

on g are exactly the structure constants C^p of g. Furthermore, the connection coefficients 
of the covariant derivative V on q Q/G associated with A take the form Ff^ = C^^AJ. 
Finally, the coefficients of the g- valued 2-form uj on M are 

The basis {cq,} for Sec(g) induces a basis {cj = ^ © 0, Cq, = © Cq,} for Sec(r). W.r.t. 
this basis, the Lie algebroid bracket (fTIJ|) is given by 

[ej,ej] = — w^-e^, [ei,eo] = Lj^e^, [^a^^p] = Cl^^e^. 

This Lie algebroid structure on r : TQ/G ~ TM©g ^ M is the so-called Atiyah algebroid 
(see also [101 El). 

Suppose now that the system is subject to some linear constraints D C TQ. In contrast 
with the previous example, there is a natural fibration of Q available, Q —>■ M = Q/G, 
so it makes sense to compare D with the vertical subspace YtTq of this fibration, rather 
than to assume that D is the horizontal distribution of a connection on some fibre bundle. 
We will follow here the approach of 7J , although there are many other [21 El IH E] ■ In U\ 
two additional assumptions were made. First, we suppose that TgQ = Dq + {VtIg)^. This 
assumption means, among others, that S = {Sq = Dg n (y7iG)q \ q G Q} is a subbundle 
of TQ, D and Vna- Further, we will also assume that D (and therefore also S) is G- 
invariant. In |6j it is proved that there exists always a G- invariant metric on Q, which we 
now assume to be fixed. Let Hg be the orthogonal complement of Sg in Dg (with respect 
to this metric), then D = H (B S , while TQ = H (B Vttg- Further, if U is the orthogonal 
complement of S in Vtcq, then TQ = H (B S (BU. Due to the above assumptions, all three 
distributions are G- invariant and thus TQ/G = H/G © S/G © U/G. Let A : TQ ^ q 
be the principal connection whose horizontal subspace at q is exactly Hg. We now use 
this connection to consider the above decomposition of TQ/G in the isomorphic bundle 
TM © g, i.e. we have 

TM®0 = aA{H/G), 

© g = aAiVTTa/G) = aA{S/G) © aA{U/G) = s © u, 
aA{D/G) = TM©s. 

Analogously to the reduction of non-constrained systems with symmetry, the Lagrange- 
d'Alembert principle which describes the equations of motion on Q, can be reduced to 
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a 'variational principle' on a reduced space (cf. 0). This new principle generates the 
reduced equations, the so-called Lagrange-d'Alembert-Poincare equations. We now show 
that these equations are Lagrange equations on a subbundle of a Lie algebroid. Using 
the above notations we now define : W = D/G ~ TM © s — M and this bundle 
is the required subbundle of the Atiyah algebroid r : V = TQ/G ~ TM © g — M. 
Choose a basis {e/} of Sec(s) and denote ej = ejca- Let (a;*,W^) be coordinates on s, 
then i : W — > V has components = 5^, i]- = 0, = and = e". Moreover, the 
components for A : W ^ TM are A* =5], Aj = 0. 

Let L G C^iy) be the reduced Lagrangian, i.e. = l{vg). If Lc is the restriction of 

L to D/G, then the Lagrange-d'Alembert-Poincare equations are 



(i dLr dLr dL (a J y, ■ de 



dtdx' dx' dv^\ dx 

with D^j = G^^e^e] and F^^ = rj„e? (see also equations (5.8.45-47) in [Tj and (4.11-4.14) 
in [7j). In jTl m [7| IH] one can find examples considered in full detail of nonholonomic 
systems with symmetry, for instance the snakeboard and the vertically rolling disk. 

5. Normal extremals in Lagrangian systems on Lie algebroids. The last example 
can be found in the theory of geometric optimal control theory. We first briefiy recall 
some general concepts from control theory. In control theory one studies dynamical 
systems that can be steered by external devices (typically representing a human input 
to the system). In dynamical systems theory, such systems are typically represented by a 
differential equation of the following type 

x{t) = f{x{t),u{t)), 

where x G R"" represents the configuration of the system and where u G R'^ represents 
the control functions. It should be clear that the dynamics of the system is completely 
determined (up to the initial condition) by the control u{t). A geometric framework for 
studying control theory is that of an anchored bundle, i.e. we assume that the configura- 
tion space is a manifold M and that the control space R'^ is the /c-dimensional fibre of a 
bundle V over M, with projection r. The analogue of the map / is a bundle map p from 
V to TM, fibred over the identity. An admissible curve is a curve v(t) in V such that 
d/ dtiriyit))) = p(v(t)). It is not difficult to see that in a local coordinate system this con- 
dition has the precise structure a control differential equation. Assume that a Lagrangian 
function L G C°°(V) is given. In optimal control theory one studies the following varia- 
tional problem: 'among all admissible curves v defined over the interval [to;^i] and such 
that r(v(to)) = Xq and r(v(ti)) = Xi for fixed endpoints Xq,Xi G M, which one minimises 
a cost functional / L{y{t))dt ?'. The maximum principle j2I] gives necessary conditions for 
admissible curves to be minimising. Locally they are given as follows: an admissible curve 
y{t) = (x*(t), v"(t)) is minimising if there is (1) a curve Pi{t) in T*M along x^{t) and a con- 
stant real number pq = 0, —1 such that {po,Pi(t), . . . ,Pn(t)) 7^ Vt and (2) the following 
Hamiltonian system is satisfied at all time t, with H{x'^, y"',Pi) = Pip'^{x^ , v") +pqL{x^, v"^) 
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a function on T*M x V: 



x\t) = — = p\x\t)y{t)), 

dH _ dpi dL 
dH dp' dL 

The latter condition says that the function t— > H{x\t),pi{t),y"-) on the fibres of V 
attains a local maximum at the point v" = v°-{t). A coordinate free version for this 
theorem was proven by H.J. Sussmann in [21|. However, for the sake of simplicity, we 
continue to work in a local coordinate system. It is very interesting to note that there 
are two kinds of solutions: those admissible curves for which there is a Pi{t) satisfying 
the conditions with pq < and those for which pq = 0. The latter are called abnormal 
extremals since in this case the conditions from the maximum principle do not depend on 
the cost function (see also O dl I2ni)- If < 0, then the solutions are called normal 
extremals. We now assume that the bundle V has the structure of a Lie algebroid, and 
that the map p is the anchor map. We will show that the above conditions from the 
maximum principle for normal admissible curves can be rewritten as the solutions to 
Lagrangian systems on Lie algebroids with constraints. For that purpose we consider the 
equation expressing the maximality condition: 

dH dT 

= —{x\t)y{t)Mt))=P^it)f^a{x'it)) - -^SxKt).y\t)). 



Since we assumed that the Lagrangian is regular, we can consider the inverse of WL = 
dL/dv"' : V V*, which gives us the following condition on the control curve: 

v«(t) = (FL-ir(x^(t),p.(t)pl(x^(t))). (21) 



If we substitute this in the Hamiltonian equations from the maximum principle, we obtain 
after some straightforward calculations (and taking into account the structure equations 
of the Lie algebroid V) that (a;*(t), v"(t)) is a solution to the Lagrangian equations on the 
Lie algebroid: 

= p^(x)v", 

lfdL\ ^.dL_dL (22) 
dt\d\/'^J '^"dx' "'^ 9^=" 

This system is in general not equivalent with the Hamiltonian system from the maximum 
principle. The condition that v° is in the image of WL~^ has to be taken into account. In 
the specific case that FL~^ is a linear map, this is precisely saying that v(t) is contained in 
a linear subbundle W, i.e. v(t) e W = im(FL-^ o i{T*M)). Solutions of the systems (EH) 
and ()22|) are in fact special solutions of the system (0)). So, normal extremals are solutions 
of a Lagrangian system on a subbundle of a Lie algebroid. However, the equations of 
motion for the normal extremals are not equivalent to the non-holonomic equations of 
motion. The normal extremals have to satisfy stronger conditions ()22|1 . 



15 



A sufficient condition for the image of FL^^ to be a linear subbundle W of V is the 
condition that L is a Lagrangian of mechanical type, i.e. when L = T — t*V, where T 
is the kinetic energy associated with a metric on V and where is a potential function 
defined on M. A typical example of a mechanical system on a Lie algebroid is the spinning 
top (cf. fTB]). 

6 Conclusions and Outlook 

In this contribution we have dealt with systems on Lie algebroids that are subject to 
some constraints. We obtained an intrinsic description of the dynamics of the systems 
in terms of a section of a prolongation bundle. Our approach unifies models for both 
Lagrange-d'Alembert equations and Lagrange-d'Alembert-Poincare equations. 

Such a simultaneous description can be very handy in many applications. For example, in 
[H] it is shown that the so-called 'Lagrange-Poincare bundles' form the ideal platform on 
which the reduction process of (unconstrained) Lagrange-Poincare equations can be re- 
peated. Knowing that such bundles are in fact Lie algebroids, we hope that Lie algebroid 
theory (and in particular the description of such systems as sections of a prolongation Lie 
algebroid) will play an important role in future developments of the process called 'La- 
grangian reduction by stages'. Due to the observations in this paper, it has become clear 
that, if we want to find a geometric formalism for successive reduction of non-holonomic 
systems, we will need to explore the geometry of (the prolongation of) Lagrange-Poincare 
subbundles. 

Also in the case of non-holonomic systems (example 3), the above developped theory leads 
to interesting new insights. Indeed, the prolongation bundle T*W —>■ W, where e.g. the 
Poincare-Cartan 2-form 66l and the dynamics P live, plays a crucial role in our approach, 
but its importance has, so far, not been recognized in the literature. A framework that 
seems to be closely related to ours is [221 12n] (although their set-up is more general since 
also time-dependent systems were included). In those papers, two fundamental two-forms 
on TW have been considered. It would be of interest to explore the relations between 
those two forms on the one hand and 66 l on the other hand. Further, we intend to find 
out how our operator 6 ffis in the approach of I23j. 

An other path for future developments is that of an appropriate framework for studying 
Hamiltonian equations on a subbundle of a Lie algebroid. Hamiltonian systems on Lie 
algebroids were already considered in e.g. fOlE]- If t* : V M and /i* : W* M are 
the duals of r and fi, then the main object of a Hamiltonian description for constraint 
systems will be a section of the prolongation bundle (fi*)^ '■ T^\N* — ^ W*. Similar as 
before, also the, not unrelated, bundles (r*)'', (t*)^ and {^*)^ will come into the picture. 
In the special example of systems with symmetry, we should be able to relate the ffist 
equation in (|2Uj) . in a Hamiltonian formulation, with the momentum equation (for a recent 
survey see [T]). 

Acknowledgements. We are indebted to Frans Cantrijn, Eduardo Martinez and Willy 
Sarlet for useful discussions. 
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